ABSTRACT 

Most of the present understanding of the S=l quantum spin chains displaying the 
Haldane gap is coming from the so-called valence-bond-solid (VBS) Hamiltonian which has 
an exactly known ground state. We show that this point is characterized by the onset of 
short-range incommensurate spin correlations in the one-parameter family of Hamiltonians 
Hq = cosOJ^i Si ■ Sj+i + sin6> ^^(S, • Si+i) 2 . This gives a physical meaning to this special 
point. We establish precise values for the gaps, correlations, the string order parameter, 
and identify the VBS point as a disorder point in the sense of classical statistical mechanics. 
It is a quantum remnant of the classical transition between a ground state with long-range 
Neel order and a ground state with incommensurate long-range order. 
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I. INTRODUCTION 

It has been conjectured by Haldane 1 that antiferromagnetic quantum spin chains have 
a disordered ground state with a gap to spin excitation when the spins are integer. This 
phenomenon has been studied extensively over the years and a simple physical picture has 
emerged through the consideration of the so-called valence-bond solid (VBS) Hamiltonian 2 . 
This peculiar Hamiltonian contains, in addition to the simplest isotropic bilinear nearest- 
neighbor exchange, a biquadratic term in the case of the spin S=l chain. This modification 
transforms the Hamiltonian in a sum of projection operators and, as a consequence, the 
ground state is known exactly and has a simple structure. This is different from the exactly 
integrable models: here nothing is exactly known about the excited states. It is believed 
that this model is smoothly connected to the usual nearest-neighbor antiferromagnet: they 
share the same physics. More precisely, there is a hidden topological long-range order 3 that 
is common to both Hamiltonians 4 ' 5 and which is revealed clearly in the VBS Hamiltonian. 
It is interesting to note that a similar situation happens also in the fractional quantum Hall 
effect 6 : here the Laughlin wave function which is the exact ground state of an approximate 
Hamiltonian does possess the hidden order which is revealed in the anyonic gauge. The 
VBS nature of the ground state of the S=l spin chain has led to the curious consequence 
of effective spins S=l/2 at the end of open chains: this has been observed theoretically 7 
by numerical means and experimentally 8 ' 9 . 

If we concentrate on the VBS model in the S = 1 case, it can be written as the sum 
of a bilinear and a biquadratic spin-spin interaction between nearest neighbors. It is thus 
natural to study it as a special case of the general bilinear-biquadratic isotropic quantum 
S = 1 chain: 

H e = cos 9 Si • S,+i + sin 9 ^(S, • S, +1 ) 2 , {I. I) 

i i 
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with 9 varying between and 2n. All energies are measured in units of the global exchange 
coupling which is omitted everywhere in this paper. The VBS Hamiltonian corresponds 
to the value 9 v b s with tan6>„b s = 1/3. In this case, each term in the sum in Eq.(I.l) 
Si ■ Sj+i + (Si ■ Si + i) 2 /3 is the projector on the spin S=2 state of the two neighboring spins 
i, This fact leads to a simple ground state wavefunction 2 . The behaviour of this model 
as a function of 9 has been studied by numerous authors 10-30 . If one increases 9 starting 
from the bilinear Hamiltonian 9 = which is known to possess a Haldane gap there is no 
phase transition till 9 = n/4 and thus the VBS Hamiltonian (9 v b s = atan(l/3) < 7r/4) is 
smoothly connected to the usual bilinear Heisenberg model. However its precise physical 
meaning has remained so far unexplained. 

In this paper, we clarify the physical meaning of the VBS point in the phase diagram of 
the family (LI) of models. Consider first the classical limit S — > oo of Eq.(I.l). For 9 = 
the ground state is antiferromagnetically long-range ordered with ordering wavevector 
q = it (Neel state) . When 9 is increased, the order becomes incommensurate when 9 > 9 C 
with tan 6> c = 1/2: the wavevector shifts from q = ix. As a consequence, the static structure 
factor S(q) has a delta peak at q = n when 9 < 9 C and a delta peak at q < n when 9 > 9 C . 
Now, in the quantum case S < oo with S integer, fluctuations wash out long-range order 
and we are left with short-range order below a characteristic correlation length £, according 
to Haldane's conjecture. The delta peak of S(q) is thus smeared and acquire a finite width 
given by 

It is important to note that, due to this finite width, the incommensurate behaviour 
cannot be seen immediately in the quantity S(q) when q shifts away from the commensurate 
position. This is best understood by considering the analytic structure of S(q) in the 
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complex q-plane. Short-range order means that singularities (poles or branch points) are 
away from the real axis at a distance ~ In the commensurate phase, the real part of 
the nearest singularity is rc and the peak of S(q) for q real is also at q = it. If we increase 
the parameter 9, at some value the real part of the leading singularity will move away 
from q = tt: this means that real-space correlations oscillate with a new period. However, 
due to the width of the peak in the structure factor, the maximum of S(q) remains at 
q = 7r till the shift Aq of the real part reaches a value Then for larger values of 

9, the structure factor will exhibit an incommensurate peak. We have obtained evidence 
that, right at the VBS point, the correlations become incommensurate in real space: it is 
a "disorder" point in the language of classical statistical mechanics 31 . For a larger value of 
9, the function S(q) exhibits the signature of incommensurability 29 : this point is properly 
called a Lifshitz point. This splitting of the classical phenomenon at 9 C is typical of systems 
with only short-range order. 

We proceed by first recalling briefly in section II the state of knowledge on the bilinear- 
biquadratic S = 1 quantum spin chain. In section III we import the concepts of so-called 
" disorder" points from classical statistical mechanics. They are discussed by use of simple 
classical spin models. In section IV, using the Density Matrix Renormalization Group 
algorithm 32 , we calculate energy gaps, correlation functions and correlation lengths, as 
well as the S = 1 string order parameter in the neighborhood of the VBS point. We 
demonstrate that the spin correlations exhibit a change of behaviour in real space right at 
the VBS point and that this point is a quantum example of a "disorder" point. Section V 
contains our conclusions. 
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II. PHASE DIAGRAM OF THE BILINEAR-BIQUADRATIC SPIN CHAIN 

Let us represent the phase diagram of models He as in figure 1. For 9 = and 9 = n, 
one finds the isotropic (anti) ferromagnetic quantum Heisenberg model. Antiferromagnetic 
(resp. ferromagnetic) models corresponds to — n/2 < 9 < +ir/2 (resp. +n/2 < 9 < 3tt/2). 
Some points in the phase diagram have been studied in detail and we summarize below 
the current knowledge: 

• = 0: Isotropic antiferromagnetic quantum Heisenberg model: this well-studied model 
has a non-degenerate disordered ground state with exponentially decaying antiferromag- 
netic correlations obeying a law (So ■ S n ) « (— ) n exp(— n/£)/y/n and a gapped spectrum 
(Haldane gap 0.41). The static structure factor S(q) is a square-root Lorentzian peaked 
at q = tv. 

• 9 = 0.1024-7T (tan# = 1/3): this is the VBS model with exact valence-bond-solid ground 
state. The spin correlations are purely exponential with a correlation length £ = (ln3) _1 
0.91. There is a gap in the spectrum (A = 0.664). 

• 9 = 0.25-7T: this is the Lai- Sutherland model, the Hamiltonian is a sum of permutation 
operators and exactly integrable by the Bethe ansatz 10 ' 11 . The ground state is unique 
and the model is critical. The corresponding conformal theory is SU(3)k=i- There are 
zero-energy modes for q = 0,±27r/3. 

• 9 = — 0.257r: The model is solvable exactly by the nested Bethe ansatz 12 ' 13 . One finds a 
critical system with a unique ground state. The conformal theory is SU (2)k=2- There are 
zero-energy modes at q = 0,tt. 
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• 9 = — 0.50/t: the physics is that of a dimerized state; the order parameter is given by the 
coefficient c 2 in the singlet-singlet correlation: 

((SA+i)(SA+i)> - ci + (-l) l - J c 2 , (II.l) 

for |z — j| — > oo. The ground state is twice degenerate in the thermodynamic limit and 
the spectrum is gapped (A = 0.17). The correlation length is given as £ = 42.2: these are 
exact results 18-21 . 

• 9 = — 0.75-7T: possible location of a continuous phase transition from a ferromagnetic to 
a dimerized phase 23 ' 25 ' 26 . 

• 9 = it: This is the isotropic ferromagnetic Heisenberg model. There is ferromagnetic order 
with gap less excitations. The ground state is the ferromagnetic state for n/2 < 9 < 57r/4. 

With these points, one constructs the following phase diagram 14 ' 17 : Starting at 9 = tt, 
one finds an ordered ferromagnetic state without gap. The ferromagnetic phase terminates 
at 9 = — 0.757T. A continuous phase transition leads to a dimerized state. A prediction by 
Chubukov 22 of a non-dimerized nematic phase seems refuted by Fath and Solyom 25 ' 26 . In 
the dimerized phase, the ground state is a singlet with a double degeneracy due to a Z 2 
symmetry breaking. The order parameter is given by c 2 in the correlation function (II.l). 
A continuous phase transition at 9 = —0.2bn leads to a Haldane phase, with a unique 
disordered ground state, exponentially decaying correlations and a gapped spectrum. This 
gapped phase ends at the Lai-Sutherland point 9 = +0.257T where a continuous transition 
leads to a phase which is possibly trimerized (see refs. 25,26 for a detailed discussion). 
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One is back to ferromagnetic phase for 9 = +0.5n. This phase diagram is displayed in 
figure 1. Up to now, the VBS point appears to be generic in the Haldane phase. 

Recently, Bursill, Xiang and Gehring 29 considered, using the DMRG, the Fourier 
transform of the spin-spin correlations, i.e. the static structure factor: 

S(q) = J2e iqn (Sn- So), (JJ.2) 

n 

In the Haldane phase, at the isotropic point 9 = 0, S(q) is a square-root Lorentzian with a 
peak at q = tt. Since parity is unbroken in the phases we discuss, we restrict the momenta 
to the interval (0,7r). It was found that the peak of the Fourier transform S(q) starts 
to move away from q = n towards q = ±2tt/3. This happens at tan# = 0.43806(4), or 
9 = 0.1314-7I-. They found also that when 9 — > 0.25-7T, the peak reaches 2n/3 in agreement 
with the period-3 zero modes that are seen at the Lai-Sutherland point. Their conclusion 
is then that there are three regions between 9 = and 9 = 0.50n: 

• < 9 < 9: short-ranged antiferromagnetic correlations. 

• 9 < 9 < 0.25"7r): short-ranged spiral order; the peak of the Fourier transform shifts 
from q = tt to q = 27r/3; the spectrum is still gapped. 

• 0.257T < 9 < 0.50/t: possible trimerized phase beyond the Lai-Sutherland phase 
transition. 

If one considers the classical limit of model (LI), there is a related phenomenon. For 
9 smaller than 9 C = arctan(l/2) = 0.1487T, the ground state is the usual commensurate 
Neel order with wavevector q = tt, while beyond this value of 9 the ground state becomes 
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an incommensurate spiral characterized by wavevector q such that cosq = — |cotan6>. 
Of course the classical ground states have long-range ordering and when going to finite 
spin values this order becomes short range. We will show in section IV that this short- 
range order naturally splits the commensurate-incommensurate transition in two distinct 
phenomena: one happens at 9 v b s where the spin oscillations becomes incommensurate 
in real space and one happens at 9 where incommensurability becomes obvious in the 
structure factor. 

Before discussing our results, we now recall the corresponding concepts of classical 
statistical physics, first developed by Stephenson 31 . 

III. SHORT-RANGE ORDER AND "DISORDER" POINTS 

If one starts from a classical model and considers finite integer spins then, according 
to Haldane's conjecture, there is only short-range order and a finite correlation length. 
This is, roughly speaking, an example of " quantum paramagnetism" . Finite spin is in a 
sense equivalent to a finite temperature. In Haldane's mapping 1 onto a nonlinear sigma 
model, the coupling constant is equal to the inverse of the spin while in the nonlinear 
sigma model describing classical two-dimensional systems the coupling is the temperature 
itself. This means that an integer spin chain has a physics which is related to that of a 
two-dimensional spin system at nonzero temperature. Since there is no long-range ordering 
in such a two-dimensional system according to the Mermin- Wagner theorem, the ground 
state of the spin chain is short-range ordered. We consider thus classical systems in their 
paramagnetic phase to understand the physics of finite-spin chains. Strictly speaking, one 
should consider classical two-dimensional systems but in fact, for our purposes, the physics 
is absolutely similar to that of three-dimensional systems above the critical temperature. 
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Let us consider a magnetic Hamiltonian that exhibits two ordered low-temperature 
phases, one with commensurate correlations and the other with incommensurate correla- 
tions. One may think for example of a square lattice of classical spins with nearest-neighbor 
exchange J\ and third- nearest-neighbor J 3 : when J3/J1 > 1/8 one destabilizes the Neel 
order and obtains an incommensurate spiral whose pitch evolves continuously. We note P 
any parameter that controls the zero-temperature phase transition (e.g. anisotropy, pres- 
sure, ratio of exchange couplings, etc.). A generic phase diagram is given in figure 2. We 
consider the case where these low-temperature phases are separated from the disordered 
paramagnetic high-temperature phase by continuous transitions (If the classical system is 
2D then the T c is zero and the reasoning is unchanged). It is clear that the short-ranged 
correlations in the disordered phase will be of variable nature: "close" to the commensu- 
rate phase, they will be commensurate; "close" to the incommensurate phase, they will be 
incommensurate. One can guess that there will be a line in this phase diagram, where the 
correlations change their behaviour; this change will be linked to correlations of very short 
range, thus to a state with a minimum of short range order. Hence the name of disorder 
line. If one moves along path (A) in the paramagnetic phase in figure 2, there should be 
a change in the correlations. 

If one considers the real-space spin-spin correlations along path (A) they will develop 
incommensurate oscillations at some point Ad. If one considers now the correlations in 
Fourier space S(q), one finds that the peaks of the Fourier transform still stays at the 
value for commensurate correlations, even though the real-space correlations are already 
incommensurate, due to the finite correlation length: the peak width is linked to It 
is only "closer" to the incommensurate phase that the peak will start to shift. This will 
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happen at a second point A L on path (A) in fig. 2. This is easy to understand by taking a 
simplified form for S(q): 

For convenience, we shift the momenta to set q x = 0. It is only for a < that S(q) 
has a double-peak structure. When a > is large enough, a 2 > 4£~ 2 , all the poles of 
(III. 1) are on the imaginary axis in complex-q space and the real-space correlations do not 
oscillate. But when a 2 < 4£~ 2 the poles have a real part and thus there are real-space 
oscillations. It is only when the real part of the poles is large enough that the structure 
factor itself displays a two-peak shape. This effect is entirely due to the finite correlation 
length £ i.e. short-range ordering and finite width of the peak in S(q). It is clear from the 
simple example above that it is only the analytic structure of S(q) that matters and not 
our peculiar eq.(III.l). As such, this is a general behaviour. 

The starting point of real-space oscillations is the disorder point. It extends in the 
plane of Fig. 2 in a disorder line D. The starting point for the double peak structure in S(q) 
is the Lifshitz point, extending in a line L in fig. 2. In experiments, one normally measures 
the structure factor in reciprocal space and will thus observe this line. It is also clear that 
the two lines must end in the multicritical point where the three phases meet, which is 
thus necessary for their existence. 

To avoid unnecessary generalizations, we take the results from an example treated by 
an RPA method in Ref.33. It is an Ising spin chain with a ferromagnetic J\ interaction 
between nearest neighbors and an antiferromagnetic J 2 between next-to-nearest neighbors. 
The RPA treatment shows that there are three regimes. The disorder temperature is: 

T D = J 2 /A\J 2 \ + 2\J 2 \. {III. 2) 

10 



One derives the following expressions for x large: 
T < T D : (S(0)S(x)} ~ e"^-( T ) . 

T = T D : (S(0)S(x)} ~ xe" fcoX with coshfcoa = Ji/4|J 2 |. 

T > T D : (S(0)S(x)) ~ e - x ^+( T ) cos(g(2» with g(T) ~ (T - lb) 1 / 2 . 

For £±(T) one finds that, on the commensurate side, the correlation length exhibits an 
infinite derivative at Td; it will typically be very small, but not necessarily a minimum or 
zero. The derivative on the incommensurate side is finite (see figure 3). This characterizes a 
disorder line of the first kind. There are two more special properties: (i) The susceptibility 
shows a particularly simple form at the disorder line, (ii) If one considers the correlation 
functions and compares them to an Ornstein-Zernicke correlation function (for x large) for 
a d-dimensional system: 

(S(0)S(x)} ~ e -*/€(r)/ r (d-i)/2 j ( /7/ 3 ) 

one sees that the correlation functions are those for d = 1, as expected for a chain, except 
at the disorder point: formally, they correspond to d = —1. Let us add that the incom- 
mensurate correlations are given by a wave vector q which shifts continuously from the 
commensurate value q = 0; the exponent 1/2 is however non-universal. 

In the spin chain problem with S=l, we are always in the paramagnetic phase i.e. 
we are following a path like (A) in fig. 2 when varying the parameter 9 in the model (1.1). 
We thus expect to cross the disorder point and the Lifshitz point that are the quantum 
remnants of the classical transition at 9 C . 
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IV. THE NEIGHBORHOOD OF THE VBS POINT 

For our calculations, we use the DMRG: see Ref.32 for a detailed discussion of the 
algorithm. We apply it to chains of a length L = 96 and keep M = 80 states. This 
is sufficient to find truncation errors smaller than 10~ 12 in the considered region. It is 
therefore not necessary to extrapolate results in M, as they are extremely close to the 
exact results. For the VBS point, we recover the exact results within machine precision. 
Due to the very small correlation lengths (typically smaller than 3) a length L = 96 is 
sufficient to obtain the results of the thermodynamic limit. From the DMRG viewpoint, 
this situation is ideal. There is however a problem with purely computational errors: 
the spin-spin correlations are, for a distance of 30 to 40 sites, of the order 10~ 13 or less. 
As they are obtained by summing small numbers below machine precision (in a REAL*8 
calculation) , they must be rejected. To judge the importance of this effect, we have adopted 
the following strategy: as the system under study is isotropic and disordered in the region 
where the Haldane and the trimerized phase meet, the spin-spin correlations must obey 
the relation 

(S+S-) = (S*S* + SfS]) = 2(5f 5?>. (JV.l) 

These two quantities are calculated independently; we reject correlations that show a 
deviation of more than a thousandth from this relation. As a matter of fact, we find that 
the correlations (S^S*) reach a minimum value which oscillates randomly around 10~ 14 , 
whereas (S^S~) continues to diminish regularly. We conclude that the values for the 
latter correlation are more precise; analyzing the calculation, we find that for the latter all 
the important contributions and weights show the same sign, whereas it changes for the 
former. 
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The key of our analysis is not to analyze the Fourier transform of the correlations, 
but to analyze them directly in real space. We will thus show that the VBS point, so far 
without special role in the phase diagram, is effectively a disorder point. It shows all the 
characteristics of a disorder point of the first kind, as described in section III. 

Consider the real-space correlations (figure 4) for some values of 9 between 0.10-7T and 
0.125tt. This includes the VBS point (6 vbs = 0.1024tt). The correlations for 9 < 0.1024tt are 
perfectly antiferromagnetic; this is most evident in a logarithmic plot of |(— l) n y / n(S n -So)|. 
For antiferromagnetic correlations, the curve shows no modulations. Above 9 v b s = 0.1024-7T, 
the logarithmic plots show oscillations with periods that become shorter for increasing 9, 
to end at a period of 3 for 9 — > 0.25n. Thus the VBS point is a disorder point. The 
correlations are no longer antiferromagnetic, but already incommensurate: this point is 
missed if one considers only the Fourier transform as in Ref.29. These modulations can be 
understood from the classical law for an incommensurate high-temperature phase (adapted 
to d= 2): 

(S n • So) « cos[q(9)n] =- = (-l) n cos[(tt - q)n] {IV. 2) 

V n ' V n 

The modulations should thus show a period tv/(tv — q), which is easier to see than the 
period originating directly from cos qn. To show that the correlation functions can be well 
described by (IV. 2), we have attempted a direct fit of our results. This fit is complicated 
by the fact that there are effectively three parameters to be controlled, the wave vector q, 
the correlation length £ and also a phase factor 0, from replacing n = i — j by (n — 4>) in 
the argument of the cosine. We find that the fit is extremely sensitive to the parameter 
values, which allows for a good fit. As an example, we take 9 = 0.1157T, sensibly below the 
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point given so far for the change of the correlations. We obtain the fit shown in figure 5, 
for cj) = 0.65, £ = 1.08 and it — q = 0.1987T. The wave vector q has already shifted by 20 
percent from the antiferromagnetic value q = it. 

For 9 closer to the VBS point, the fit is made more complicated by the errors in the 
correlation function due to the finite precision of the computer. To estimate the behaviour 
q{9), we therefore consider the periodicity: For a period p, n — q rc/p. The discrete nature 
of the problem limits this approach. We find the behaviour shown in figure 6. Clearly, is 
is compatible with q oc (6 — 9 v b s ) a with < a < 1, the behaviour of a disorder point of 
the first kind. The curve would be compatible with a ~ 1/2; but we are in no position to 
give a precise estimate. 

We have also calculated the Fourier transform of the correlation function and find 
results in agreement with those given by Bursill et al. 29 . The VBS point has no special 
significance for its behaviour; the point 9 = 0.1314-7T where the peak starts to shift, can 
now be identified as a Lifshitz point: see fig. 7. As expected in a system with short-range 
order, it is distinct from the disorder (VBS) point. 

For the correlation lengths, we find that they show a minimum for the VBS point, with 
an infinite slope (numerically: very large) for £ in the commensurate regime (9 < 0.10247r), 
and a slow increase in the incommensurate regime with a finite slope, given in figure 8. The 
correlation lengths have been found by different methods: in the regime 9 < 9 v b s , we have 
compared the spin-spin correlations numerically and graphically to a law exp(—n/^)/^/n, 
which were in all cases in good agreement. We estimate the precision of the results of the 
order of 1 per cent: the truncation errors are of the order 10 -13 , a serious underestimation 
can therefore be excluded. For L = 96 and £ 1 — 2, finite size effects are of no importance. 
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The situation is more complicated for the regime 6 > 6 v b s . As we have seen, the fit of the 
theoretically expected behaviour to the found curve is rather complex. If one considers 
a plot of |(— l) n y/n(Si ■ Sj)\ (figure 4), one finds that a linear fit for the maxima is quite 
good. This is stable in the sense that a factor cos qx influences the logarithm least when 
it is close to 1. Very generously estimated, the error of the graphical evaluation should 
be below 5 per cent. We estimate that for 9 > O.lbn the underestimation due to a non- 
negligible truncation error dominates. For 6 v b s < < 0.1 lir, we could not obtain the 
correlation length: on the one hand, the periods due to the incommensurability are too 
long to separate them well from the exponential behaviour. On the other hand, neither a 
fit exp— r/£ nor exp(— r/^)/y/r is satisfactory. We don't know whether this observation is 
due to the crossover of the behaviour of the correlation function or simply due to problems 
of the numerical method. In any case, the results indicate strongly an extrapolation to the 
VBS point with a finite slope. 

The other important quantities are the gap and the string order parameter O n (i,j) = 
(Sf exp(z7r J2k=]+i Sk)Sj)- The gap shows a maximum for 9 m 0.123/t, which is thus 
linked neither to the disorder nor the Lifshitz point: see figure 9. For the transition points 
8 = ±0.25/T we have not obtained serious estimates: the critical fluctuations imply a 
greater M , and the vanishing gap is difficult to see. Our results are well compatible with 
a zero gap, but not precise enough to give a serious estimate. The point 9 = — 0.20n is 
sufficiently close to the transition to cause the same problem. For the VBS point, we find 
a gap value A = 0.664, in agreement with Ref.25. Since the gap is smooth at 9 v b s , this 
implies that the spin wave velocity has a singularity at this point (c = A£). 
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The string order parameter in the thermodynamic limit i — j ' — > oo has its extremum 
for the VBS point: \O n \ = 4/9: see figure 10. One sees that the VBS point, though a point 
of minimal spin order, is a point of maximum hidden topological order. 

To complete the identification of the VBS point as a disorder point, we note that 
there is the equivalent of the particularly simple form of the susceptibility of the classical 
model of section III. The exact correlations at the VBS point obey a one dimensional 
Ornstein-Zernicke form: they are purely exponential (no prefactor), whereas the non- 
linear sigma model yields two-dimensional correlation laws. This "dimensional reduction" 
is accompanied by a particularly simple form of the Hamiltonian: it can be decomposed 
into a sum of local projection operators. The problem loses its quantum character and 
turns into a classical one-dimensional problem. 

V. CONCLUSION 

We have shown that the VBS point 2 is a disorder point in the sense of classical statis- 
tical mechanics. It is thus identified as a point which is not just by chance exactly solvable, 
but shows this property for more profound physical reasons. The results of Bursill et al. 
who have considered the Fourier transform of the spin-spin correlations fit naturally in the 
picture we have developed: the point they identify as the point where correlations change 
is simply the Lifshitz point following the definition given above. Quantum fluctuations in 
the integer spin chain wash out long-range order. As a consequence, the transition that 
happens at the classical level for 9 C is no longer a phase transition. However the change 
of short-range correlations still happens in the quantum system at the VBS point in real- 
space. Due to the finite correlation length this is not seen immediately in S(q), hence the 
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Lifshitz point which is distinct. Since £ — > oo when the spin increases one expects that 
these two points should merge in the classical limit right at 9 C . This is summarized in 
figure 11. 

This implies the following description of the phase diagram: There are three regions 
for < 6 < 0.25tt: 

< 9 < 6 v bs- Short-range antiferromagnetic correlations; Haldane phase with gapped 
spectrum. A generic description is given by the VBS model or also the isotropic antiferro- 
magnetic Heisenberg model. 

9 v bs < < 9: Incommensurate short-range correlations with a wave vector q < n, 
that shifts away from n as n — q oc (9 — 9 v b s ) a , & ~ 1/2. In the Fourier transform, the peak 
stays at q = it. The spectrum is gapped: we expect that the low-lying Haldane modes are 
now at the incommensurate wavevector. Part of the VBS physics remains valid: the gap, 
the hidden order, and the free spins 1/2 at the ends of an open chain. 

9 < 9 < 0.25-7T: The physics is similar to that of the former region; the peak of the 
Fourier transform shifts from q = it to q = ±2w/3 and the incommensurate correlations 
become visible. The spectrum is gapped. There is however no profound physical difference 
between this region and the former one. 

The above picture does not challenge conventional wisdom in the sense that the usual 
Heisenberg Hamiltonian 9 = shares the same physics with the VBS Hamiltonian and 
there are no additional phase transitions between the point 12 ' 13 9 = — 0.25-7T and the Lai- 
Sutherland point 9 = +0.257T. However the VBS point itself means that the physics has 
changed beyond 9 v b s : due to the incommensurate correlations, it is no longer possible 
to capture the low-lying excitations by a nonlinear sigma model with symmetry breaking 
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pattern 0(3)/0(2) as used originally by Haldane. In the incommensurate regime, the full 
rotation group is broken down. Appropriate nonlinear sigma models have been contem- 
plated before 34 . Since they involve a nonabelian symmetry, they are generically massive, 
consistent with the gapped nature of the Haldane phase. 

There is an interesting relationship 6 with the fractional Quantum Hall effect. The 
present understanding 35 of the physics of the FQHE at filling v = 1/m, m odd, is based on 
Laughlin's wavefunction ifj m which is the exact ground state of a truncated Hamiltonian. In 
addition, this function embodies the hidden long-range order that is revealed in the anyonic 
gauge. This is similar to the situation of the VBS wavefunction. Expectation values 
computed with the Laughlin wavefunction correspond to a classical statistical problem 
which is the two-dimensional one-component plasma (2D OCP): this is Laughlin's plasma 
analogy. The case of the full Landau level ipi corresponds to the special point 36 V = 2 
of the 2D OCP (T being the ratio of the squared electric charge to the temperature) at 
which the density correlations begin to oscillate in real space, a precursor phenomenon of 
the crystallization that occurs at V 140 when the plasma is dilute enough. This special 
point has also some of properties expected from a disorder point 36 . This is similar to what 
we observe at the VBS point. 
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FIGURE CAPTIONS 



Figure 1: 

Phase diagram of the bilinear-biquadratic S = 1 isotropic quantum spin chain as a function 
of 9. Solid lines: transition points; dashed lines: other special points. H(AFM): isotropic 
antiferromagnet; VBS: the valence-bond-solid model; L: the crossover point studied in 
Ref.29; H (FM): isotropic ferromagnet. 

Figure 2: 

Schematic phase diagram: a disordered high-T phase is linked by two continuous tran- 
sitions to two ordered low-T phases. P is a parameter which controls the nature of the 
ground state. The dashed lines represent the disorder (D) and Lifshitz (L) lines, where 
the behaviour of the correlations changes in real and in Fourier space respectively. 

Figure 3: 

Correlation length at a disorder point of the first kind (schematical from RPA): the com- 
mensurate phase is to the left. 

Figure 4: 

Real-space spin-spin correlations as a function of the distance, for several values of below 
and above the VBS point (at 6 v h s = 0.1024/t). The modulations appear above the VBS 
point. Here K n = ln[|(-)V^<S • S n )|]. 
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Figure 5: 

Comparison between the spin-spin correlations predicted (solid line) and calculated nu- 
merically (squares + dashed line) for 9 = 0.1157T, above 9 v b a , but below 9. The dotted line 
is {—) n y/nex^{n/^){SQS^) and the solid line is cos [(n — <j))(ir — q)]. 

Figure 6: 

Wave vector it — q characteristic of the spin correlations. There is a singularity at the VBS 
point consistent with the identification of a disorder point. 

Figure 7: 

The Fourier transform S(q) as a function of momentum for various values of 9. The 
two-peak structure appears only for 9 > 9 = 0.13147T while nothing is seen at 9 v b s - 

Figure 8: 

Correlation lengths for various values 9. The minimum is at the VBS point. 
Figure 9: 

Gaps in the Haldane phase as a function of 9. The maximum value is between the VBS 
point and the Lifshitz point. 

Figure 10: 

String order parameter \O n \ in the thermodynamic limit in the Haldane phase. The max- 
imum is exactly at the VBS point. 

Figure 11: 

Schematic phase diagram of integer spin chains with biquadratic coupling: model (1.1). 
Here the inverse spin plays the role of a temperature: the S=l case corresponds to short- 
range ordered phases. Of course this picture may be altered by other types of ordering 
such as dimerization. 
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